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From a macroscopic theory of the quantum vacuum in terms of conserved relativistic charges (generically 
denoted by q^""' with label a), we have obtained, in the low-energy limit, a particular type of f{R) model 
relevant to cosmology. The macroscopic quantum- vacuum theory allows us to distinguish between different 
phenomenological f{R) models on physical grounds. 
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1. INTRODUCTION 



The development of high-energy physics and cosmol- 
ogy over the last years has led to the realization that, 
most hkely, Einstein's theory of gravity needs to be 
modified. At the moment, we are only able to calculate 
higher-order curvature corrections to the Einstein action 
coming from scalar, spinor, and vector quantum fields 
propagating over a fixed classical background (see, e.g., 
Refs. [1, 2]). The full theory of the modificd-Einstein ac- 
tion requires the underlying microscopic theory - what 
is usually called quantum gravity theory. Since the latter 
theory is not yet established, phenomenological models 
are needed. Among these are the so-called f{R) models 
(see, e.g., Refs. [2, 3] for two recent reviews) which have 
ad hoc powers of the curvature invariants added to the 
linear term of the Einstein-Hilbert action [4] . 

We have proposed another approach to modified 
Einstein gravity, which is based on the treatment of 
the Lorentz-invariant quantum vacuum as an extended 
self-sustained system characterized by a conserved rela- 
tivistic charge q [5, 6]. Here, g is a microscopic variable 
describing the physics of the deep (ultraviolet) vacuum, 
but its thermodynamics and dynamics are described by 
macroscopic equations, because g is a conserved quan- 
tity. This quantity q is similar to the particle density 
in liquids, which describes a microscopic quantity - the 
density of atoms - but obeys the macroscopic equations 
of hydrodynamics, because of particle- number conser- 
vation. Different from known liquids, the quantum vac- 
uum is Lorentz invariant. The quantity q must, there- 
fore, be Lorentz invariant. This treatment has allowed 
us to discuss both the thermodynamics and the dynam- 
ics of a Lorentz-invariant quantum vacuum. 
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The thermodynamic approach [5] assumes that the 
quantum vacuum is a stable self-sustained equilibrium 
state, which is described by compressibility and other 
characteristics of the response to external perturbations. 
In this approach, we have found that the vacuum energy 
density appears in two forms. 

First, there is the microscopic vacuum energy den- 
sity which is characterized by an ultraviolet energy 
scale Euy, so that e{q) ^ ^ijv Most hkely, Euy 
corresponds to the standard Planck energy scale de- 
fined in terms of Newton's gravitational constant Gn, 
^Planck = ^/hcFJG^ « 1.2 X 10^** eV. But, here, we pre- 
fer to stay as general as possible and keep Euy distinct 

from £^pianck- 

Second, there is the macroscopic vacuum energy 
density which is determined by a particular thermo- 
dynamic quantity, evac(9) = e — qde/dq^ and it is this 
type of energy density which contributes to the effec- 
tive gravitational field equations at low energies. For 
a self-sustained vacuum in full thermodynamic equilib- 
rium and in the absence of matter, the effective (coarse- 
grained) vacuum energy density Cvac (?) is automatically 
nullified (without fine tuning) by the spontaneous ad- 
justment of the vacuum variable q to its equilibrium 
value qo, so that evac('?o) = 0. This implies that the 
effective cosmological constant A of a perfect quantum 
vacuum is strictly zero, which is consistent with the re- 
quirement of Lorentz invariance for zero external pres- 
sure. 

The dynamic approach [6] demonstrates how, in 
a fiat Friedmann-Robertson-Walker universe, the vac- 
uum energy density ?vac (effective cosmological "con- 
stant") relaxes from its natural Planck-scale value at 
early Planckian times to a naturally small value at late 
times. 

In this Letter, we show that the macroscopic theory 
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of the quantum vacuum, when applied to cosmology, 
gives rise to a specific class of f{R) models of modified 
gravity. 



2. GRAVITY WITH F-FIELDS AND 
LOW-ENERGY MATTER 

We consider the general case of several conserved 
microscopic variables q^'^\ for a = 1, ... , n, and 
corresponding chemical potentials /i^") [7, 8]. As in 
Refs. [5, 6], each variable q'-°'> can be represented by 
a four-form field i^("): 

(p(a)\2 ^ p{a) p(a)fiupiy p(a) _ rrj 4(0.) 

\ ; — 24 til^pcr ' /Jl^pcr — ^ [fj.^'-vpcr] " 

(1) 

The action of the four- form fields F^'^'>{x), the matter 
field (t){x), and the gravitational field g^^(x) is given by 



SiA^'^Kg,^] 



?i<^) 



X [KiF^'^)) R + e{F^''>,^) + -d^cl)d^'<j) 



j4 / 

a X \i —g 



(a) 



(2) 



For simplicity, low-energy matter is represented by a 
single real scalar field (/). The generalized potential 
e{F^°'\(j)) includes self-interactions and interactions be- 
tween all fields F^'^^ and (j), but contains no derivatives 
of these fields and explicit factors of the metric field 
5^^ or its inverse. The gravitational coupling parame- 
ter K is determined by ultraviolet physics and, there- 
fore, depends on the microscopic vacuum variables i^'°'. 
Here, and in the following, we adopt the conventions 
of Ref. [4], in particular, those for the Riemann tensor 
and the metric signature ( — h H — h). Natural units with 
h = c = 1 are used throughout. 

The generalized Maxwell and Klein-Gordon equa- 
tions from action (2) read 



V„ 



p(a)p.vpa / Q^ 



F(«) 



9F(°) 



R 



dK 



= 0, (3a) 



de 



n(/) - — - == , (3b) 



where D denotes the invariant d'Alembertian operator 
and the partial derivatives d/dF'-"-^ and d/dcf) stand for 
pointwise differentiation. The variation of (2) over the 
metric g^^ gives the generalized Einstein equations: 



2K (r^, - i Rg^^ + Rg^, ^ F^'^ 



-2(\/f,\/,-g^,ajK-€{F 



a=l 



dK 



.9p. + C = 0' (4) 



with effective potential 



^^(.) 



b=l 



de 

dFib) 



(5) 



and scalar-field energy-momentum tensor 

1 



T. 



j:i^d^^d,4^-^g^,dp<f>df<f>. (6) 

Using F^'^^^'^P" as given by (1), the Maxwell equa- 
tions (3a) can be rewritten in the following form: 



da 



de 



R 



dK 



^ \dF('') ' ^^ dFi'^ 
The solution of these 4n equations is 



de 



i^^^^.H 



(7) 



(8) 



where the /i'"' are n integration constants. Eliminating 
dK/dF^"^^ from (4) and (8), one finds for the generalized 
Einstein equations 



-2K[R^, - \Rgpu) - 2(V^V, - g^, Dj K 

n 

+ (e-Y^^,^^^F^'^))g,^^T^^. (9) 



Equations (8) and (9) can also be obtained if we use, 
instead of the original action, the following effective ac- 
tion: 



^eff[^("\M^°\.9,' 



d^x^r 



9 



X KR + e 



n \ 

a=l / 



(10) 



The ^(")f("' terms in this action do not contribute to 
the equations of motion (3a), because they are total 
derivatives, 

/ d'x VWlf^^"' ^^"^ = -^ e^''" j d'^ Fi%^ . 

(11) 

The constant y^^"-* is seen to play the role of a Lagrange 
multiplier related to the conservation of the vacuum 
"charge" F^"^). 

Instead of the large microscopic energy density 
e{F^ \ 4>) in the original action (2), a potentially smaller 
macroscopic vacuum energy density enters the effective 
action (10), namely. 



pv = e{F 



(a) 



^^(a)^(a) 



(12) 



o=l 



Precisely this macroscopic vacuum energy density grav- 
itates and determines the cosmological term in the grav- 
itational field equations (9). 
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3. EQUILIBRIUM VACUA AND STABILITY 
CONDITIONS 

The main goal of our approach is to describe the 
thermodynamics of the cquihbrium vacuum [5] and to 
consider cosmology as the dynamics of relaxation to- 
wards an equilibrium state [6, 9, 10]. That is why, in 
what follows, we assume that the universe is close to 
equilibrium and that all its parameters, including the 
fields i^("^ and the chemical potentials /i^"-*, are close to 
their equilibrium values. A static homogeneous equilib- 
rium vacuum, in the absence of thermal matter, corre- 
sponds to a stationary point (i^o , "^o) ofEqs. (3b), (8), 
and (9) for i?^, = T^ = 0: 



de 



de 



dF(^) 



M) 



E 

a=l 



^(a)^(a) ^Q^ (^3) 



where the last equation demonstrates that the vacuum 
energy (12) is zero in an equilibrium vacuum, pvloq — 0. 

One can see the difference between the conventional 
matter field (f) and the conserved vacuum fields F^""' , 
as only the fields F'-°-'> provide the integration constants 
^w which arise dynamically from the solution (8) of the 
generalized Maxwell equations (7). (The field equations 
of generic matter fields do not give rise to such integra- 
tion constants.) These integration constants /x^"-* play 
the role of chemical potentials in thermodynamics and 
are thermodynamically conjugate to the density of the 
conserved quantities F^'^' . With appropriate nonzero 
chemical potentials, the large vacuum energy e{F^°'\ (f) 
is reduced to pv = in a homogeneous equilibrium 
vacuum state according to (13). Specifically, two large 
quantities, e{F'-°'\(j)) and X]a=i A'''"''-P''"°\ each of order 
Efjy, cancel each other due to the self-tuning mecha- 
nism [5]. This is the main property of a self-sustained 
vacuum. 

For the case of a single vacuum variable F, the chem- 
ical potential fi in equilibrium is completely fixed by the 
constraint py = 0. But, for the case of several vari- 
ables F'^°'\ there are n — 1 degrees of freedom, since the 
equation py = gives only a single constraint on the n 
chemical potentials /i'"^. This allows for the existence of 
many different equilibrium vacua and also for the coex- 
istence of several vacua [7, 8] . This last observation may 
give microscopic support to the multiple point principle 
which postulates the existence of a number of phases 
with the same energy density (see, e.g., Ref. [11] and 
references therein). 

The stationary point {Fq ,(j)o) of the thermody- 
namic potential W = e{F'-''\(j)) - ^^ ^('■)f('^) [i.e., 
the solution of (13)] should correspond to a minimum. 



which can be local or global. In particular, the vacuum 
compressibility introduced in Ref. [5] must be positive: 



Xo 



V" p{a)p{b) 



a, 6=1 



gp{a)Qp{b) 



>0. 



pial^pM 



T',4>=<t>o 



(14) 

Furthermore, the effective Newton's constant must be 
positive in an equilibrium vacuum: 



G^ 



1 



167ri^(F(|"^; 



>0, 



(15) 



in order to have a physically consistent description of 
an attractive gravitational force. More specifically, neg- 
ative K gives the wrong sign of the kinetic term for the 
graviton (which becomes ghostlike) and the quantum 
vacuum is unstable. As g-theory is based on a stable 
self-sustaining vacuum, having K > is a necessary 
condition for stability of the vacuum, together with the 
compressibility condition (14). 



4. f{R) MODEL FROM q-THEORY 

Modified-gravity f{R) models appeared already in 
the 1960s (see, e.g., Refs. [12, 13]) and were used to con- 
struct an inflationary model of the early universe in the 
1980s [14]. More recently, these f{R) models have re- 
ceived renewed attention as a way to explain the inferred 
cosmic "dark energy" by attributing it to a modification 
of Einstein gravity (see, e.g., Refs. [2, 3, 15, 16] and ref- 
erences therein). These models are, in fact, purely phe- 
nomenological models, which, in their simplest form, re- 
place the linear function of the Ricci scalar R from the 
Einstein-Hilbert action term by a more general function 
f{R)- This function /(i?) can, in principle, be adjusted 
to fit the astronomical observations and to produce a 
viable cosmological model. 

To obtain f{R) from g-theory [6], one can express 
_F(°^ in terms of R, (j), and /z*^"-* by use of the equa- 
tion system (8) and substitute the resulting functions 
F("'(i?) into (9). This can be done in a general way 
(see Ref. [3]); but, since we consider the relaxation to an 
equilibrium vacuum, we are only interested in the sim- 
pler situation of a system already close to equilibrium. 
In addition, we will restrict ourselves to a single i^-field 
(the generalization to n fields is straightforward) and we 
also omit the explicit matter (/(-field, keeping only the 
general matter energy- momentum tensor T^. 



For a single F-field, (8) gives 



dejF) 
dF 



-R 



dK{F) 
OF 



(16) 
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Close to the equilibrium state determined by (13), one 
can expand the microscopic variables as follows: 



F ^ Fo + 5F, fi^ fiQ + Sfi. 



(17) 



Expressing SF in terms of R and 6)1 and excluding 6F 
from the Einstein equations (9), one obtains: 



-2K[R^, - \Rg^,,.) + 2x (V^V, R - g^, U R 
in terms of the new dimensionless parameter 

\ / F=Fo 



(18) 



(19) 



In (18), we have omitted the ((5/i)^ term and kept only 
the leading term containing 5^. Expanding K in the 
first term of (18), K ^ (IBttGn)"^ + 5F dK/dF, one 
obtains the following modified Einstein equations: 



Rf,u - ^Rg^^u) + 2 x( V^V^ i? - 5^^ D i? 



87rG'N 

T R g^iv + R Rfiiy] - Fo 5^ g^v = T^^, , 



(20) 



where we have kept only the leading term Fq 5^ and 
omitted terms R5^ and (5/i^. 

The field equations (20) correspond to the following 
phenomenological model: 



•^phcnom — / Ci X 



1 



IGttGn 
f{R) = -R+ IBttGn X i?^ - 2 a , 



f{R) + Z:^ ) , (21a) 
(21b) 



with C^^ the standard matter Lagrange density. The 
function f{R) found belongs to the class of models 
f{R) ^ —R + i?^/(6M^), where a bare cosmologi- 
cal constant may or may not be added and M = 
{ZcPf/dR'^)^^/'^ is the scalaron mass (see Refs. [14, 15] 
and references therein; note that our sign convention for 
R is opposite to that in Ref. [15]). From q theory, the 
scalaron mass square is given by 



M2 = l/(967rGNx) 



(22) 



with X defined by (19) for a single "charge" F and 
by (25) below for multiple "charges" F^"^. At this 
moment, we have two parenthetical remarks. First, 
the effective action for (20) does not contain Rf^i.R^" 
or RftupaR^'^'"^ terms in addition to the R^ term of 
(21b). Second, a purely phenomenological connection 
between generalized-cquation-of-state models and f{R) 



modified-gravity models has also been noted in App. A 
of Ref. [17]. 

In g-theory, the cosmological constant A in (21b) 
is induced by the deviation of the chemical potential 
fi from its equilibrium value (iq: A = 87rGN pv = 
— 87rGN-Fo(5/i. For the case of n charges F^'^\ the cos- 
mological constant A is given by 



A = 87rGN PV = -87rGN Yl -^o '^M 



^''hni-) 



(23) 



For the general case of a quantum vacuum described by 
conserved vacuum variables g*-"-*, one has: 



A 



A 



87rGN 



= Pv 



-^,(«)5^ 



(a) 



(24) 



This last equation also follows from the zero- 
temperature Gibbs-Duhem relation [5] applied to 
the thermodynamic system characterized by several 
conserved variables: /x'"-* is the variable thermodynam- 
ically conjugate to the conserved variable q'-°-\ 

Note that f{R) phenomenology adds a pair of ther- 
modynamically conjugate variables, R and K . This fol- 
lows from, e.g., the thermodynamic potential e + K R — 
^^(a)^(a) jj^ (10). The corresponding thermodynamic 
identities can be used to obtain, for example, the dimen- 
sionless parameter x in (21b), which can be interpreted 
as —dK/dR evaluated at i? = 0. One then has, at 
i? = 0. 



E 

a,b 

E 

a.b 



'dR ^ ^ ^ aF(") dR 

a 

dK dK dF^^^ 



E 



dK dK 



9F(°) 9^(") 



Qp(a) Qp(b) Q^(a) 



dK dK 



dFi") dFC') \dF('')dF('>'> 



(25) 



which gives, using definition (14), precisely (19) for the 
case of a single charge F. However, this thermodynamic 
description is only applicable to f{R) phenomenology, 
since it does not take into account non-/(i?) terms in the 
action such as i^^^i?^"^ and Rfj.i.paR'^'^'"^ which appear 
due to quantum corrections [1, 2]. 

For g-theory of the type considered in (2), there are 
two stability conditions, namely, condition (15) for the 
effective Newton's constant Gn and condition (14) for 
the zero-temperature vacuum compressibility xo- These 
two conditions of g-theory correspond, respectively, to 
the following stability conditions of f{R) models (see, 
e.g., Eq. (7) of Ref. [15]): 



-fiR) >0, fiR) 



R=0 



>0. 



R=0 



(26) 
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where the prime denotes differentiation with respect to 
R and the extra minus sign traces back to our conven- 
tions, as mentioned aheady a few hues under (21b). 
These conditions ensure the stabihty of the empty flat 
universe, which is a basic assumption for g-theory as 
it concerns the thermodynamics of the equihbrium vac- 
uum. In f{R) models, the second condition of (26), 
i.e., the positive mass squared of the scalaron, M^ > 0, 
implies the stability of the Minkowski vacuum, whereas 
for the negative mass squared, M^ < 0, the Minkowski 
vacuum experiences the scalaron instability [15]. 

For the case of M^ > 0, the relaxation of the uni- 
verse to equilibrium is accompanied by oscillations of R 
with frequency M (see Ref. [15] for the f{R) model and 
Ref. [6] for the q-theory which also has SF oscillations). 
In particular, the vacuum energy density in q-theory 
has been found to relax as follows [6]: 



Pv 






sin^ Mt . 



(27) 



for cosmic time t 3> l/M. While the q field itself has an 
equation of state parameter w = — 1, corresponding to 
vacuum energy density and a cosmological constant, the 
decay of the vacuum energy density in (27) simulates the 
evolution of cold dark matter with w = 0. The effective 
equation of state w = is induced by the interaction of 
the q field with gravity [6]. 



As a result, the dimensionless parameter x in (19) can 
be expected to be of order unity. If the assumption of 
a single fundamental energy scale holds true, one then 
obtains that the mass parameter (22) of the induced 
f{R) = -R + i?V(6M2) model is of the order of the 
ultraviolet scale, M ^ Etjy- In this case, the R^ cor- 
rection is of the order of the quantum R^ corrections to 
the Einstein action [1, 2]. Moreover, as g-theory starts 
from the assumption that the vacuum is a stable self- 
sustained medium, one has M^ > 0, as follows explicitly 
from (14), (19), and (21b). 

All this agrees with the simplest ,f(R) ^ —R + 
R^ /M"^ model having a Planck-scale mass M [14], but 
disagrees with certain other models which use more 
complicated phenonienological functions f{R). For the 
particular models suggested in Ref. [15], the functions 
f{R) expanded around R = have negative M^ (mak- 
ing Minkowski spacetime unstable) and — M^ of the or- 
der of A now, where A now = SttGn pv.now ~ (10~^^ eV)^ 
is the present (positive) cosmological constant as de- 
termined by observational cosmology [18]. These cos- 
mologically desirable f{R) models can, in principle, 
be obtained by using special choices of e(F*^°\ </)) and 
K{F'^°-^). But such choices require a careful tuning of 
the parameters, which is unnatural from our point-of- 
view on the properties of the self-sustained quantum 
vacuum. 



5. COMPARISON WITH 
PHENOMENOLOGICAL f{R) MODELS 

As explained above and, more briefiy, in Ref. [6], F 
theory (or, more generally, g-theory) may give a mi- 
croscopic justification for the phenonienological f{R) 
models used in theoretical cosmology and may allow 
for a choice between different classes of model functions 
f{R) based on fundamental physics. Close to equi- 
librium, the effective ,f{R) model emerging from the 
g-theory of quantum vacuum belongs to the class of 
f{R) ^ -R + R^/M"^ models. This rules out so-called 

l/R models, i.e., models with f{R) R + M'^/R (see, 

e.g., Refs. [2, 3, 16] and references therein). 

Here, we have assumed that the physics of the vac- 
uum variable in g-theory [5, 6] is determined by a unique 
microscopic energy scale Euy. This implies the follow- 
ing orders of magnitude: 



Ko 1 Ipo 1 


■" Euy , 


(28a) 


f^dK/dFl,''^\^G^'r 


" -E-uv > 


(28b) 


HfI,''\M\^ 


" -E-uv > 


(28c) 


Xo - 


■"^uv- 


(28d) 



6. CONCLUSION 

In g-theory [5, 6] with chemical potentials /i'"' at 
their equilibrium values Hq , the vacuum energy den- 
sity pv has been found to relax, according to (27), from 
its natural Planck-scale value at early times when the 
system is far from equilibrium to a naturally small value 
at late times when the system is close to equilibrium. 
This solves the main cosmological problem: the present 
cosmological constant is small compared to Planck-scale 
values simply because the universe happens to be old 
compared to Planck-scale times. The remaining prob- 
lem is to understand why the cosmological constant does 
not completely relax to zero as i ^ oo or, in other words, 
to determine the origin of the small residual part of the 
vacuum energy which remains (almost) constant during 
the present epoch. 

The suggested f{R) model in which the cosmolog- 
ical constant appears just at the latest stage [15] does 
not follow from the macroscopic quantum-vacuum ap- 
proach. That is why one needs to find another explana- 
tion for the observed value of A. In order to produce a 
nonzero A, one must find, according to (23)-(24), pro- 
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cesses which shift the chemical potentials /x*^"-* away from 
the equilibrium mean-field values H(^ ■ 

There has been the suggestion to relate the present 
value of A to quantum fluctuations of the vacuum en- 
ergy, for the case that the vacuum energy itself is nulli- 
fied (see, e.g., Ref. [19]). This approach could be helpful 
in g-theory, where the vacuum energy is necessarily zero 
in full equilibrium. The observed nonzero value of A 
would then correspond to quantum or thermodynamic 
fluctuations of the chemical potentials y^*^"-* compared to 
the equilibrium mean-field values /Xq . Still, there could 
be other processes which are able to shift fi^"-^ somewhat 
away from the mean-field values /Iq , one example be- 
ing the backreaction of matter and another example be- 
ing matter production by the rapid oscillations which 
accompany the relaxation of the vacuum energy (27). 
It remains to be seen which of these processes, if any, 
is the dominant one for the inferred small but nonzero 
value of the vacuum energy density from observational 
cosmology. 
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